Pion form factor from a contact interaction 
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In a Poincare-covariant vector-boson-exchange theory, the pion possesses components of pseu- 
dovector origin, which materially influence its observable properties. For a range of such quantities, 
we explore the consequences of a momentum-independent interaction, regularised in a symmetry- 
preserving manner. The contact interaction, whilst capable of describing pion static properties, 
produces a form factor whose: evolution for Q 2 > 0.17 GeV 2 disagrees markedly with experiment; 
and asymptotic power-law behaviour conflicts strongly with perturbative-QCD. 
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The pion has a unique place in the Standard Model. It 
is a bound-state of a dressed-quark and -antiquark, and 
also that almost-massless collective excitation which is 
the Goldstone mode arising from the dynamical break- 
ing of chiral symmetry. This dichotomy can only be un- 
derstood by merging the study of many-body aspects of 
the QCD vacuum with the symmetry-preserving analysis 
of two-body bound-states. Furthermore, the possibility 
that this dichotomous nature could have wide-ranging 
effects on pion properties has made the empirical inves- 
tigation of these properties highly desirable, despite the 
difficulty in preparing a system that can act as a pion 
target and the concomitant complexities in the interpre- 
tation of the experiments; e.g., [1-4]. 

A true understanding of the pion is achievable via 
QCD's Dyson- Schwinger equations (DSEs) [5]. Key to 
this is the existence of a nonperturbative symmetry- 
preserving truncation scheme [6-8], which enables: the 
connection to be made between dynamical chiral symme- 
try breaking (DCSB) and the bound-state pion [9]; and 
the use of experiment to chart QCD's /3-function. The (3- 
function's perturbative evolution is well-known but this 
new feedback between experiment and nonperturbative 
theory may provide access to its long-range behaviour. 

Studying the electromagnetic pion form factor, 
F° ni (Q 2 ), is an ideal way to elucidate the potential of 
this interaction. For example, there is a prediction [10- 
12] that Q 2 F-„{Q 2 ) « constant for Q 2 ^> m 2 in a theory 
whose interaction is mediated by massless vector-bosons. 
The verification of this prediction is a strong motivation 
for modern experiment [1-3]. However, one may adopt a 
different view of this programme; namely, as an attempt 
to constrain and map experimentally the nature of the 
exchange interaction that binds the pion. 

Poincare covariance entails that the Bethe-Salpeter 
amplitude for an isovector pseudoscalar bound-state of 
a dressed-quark and -antiquark takes the form 

r£(A; P) = r J 75 [iE n (k; P) + 7 ■ PF n {k; P) 

+ 1 -kG^{k;P) + a llv k ll P v H v {k;P)] , (1) 

where k is the relative and P the total momentum of the 



constituents, and {r J ',j = 1, 2, 3} arc the Pauli matrices. 
(We employ a Euclidean metric with: {7^,7^} = 25 



it 



7m! 75 = 74717273; and a - b = I]i=i a «M Tnis 
amplitude is determined from the homogeneous Bethc- 
Salpeter equation (BSE): 



n(k;P)} 



(27T)" 



■ixl(q;P)]srKZ:( q ,k;P), (2) 



where P) = S(q + P)rj(g; P)S(q), r,s,t,u rep- 

resent colour, flavour and spinor indices, and K is the 
quark-antiquark scattering kernel. In Eq. (2), S is the 
dressed-quark propagator; viz., the solution of the gap 
equation: 

S(p)~ 1 = «7 • p + m 



d 4 q 



g 2 D^(p- q)—j IJ/ S(q)—r i/ (q,p), (3) 



wherein m is the Lagrangian current-quark mass, D^ v (k) 
is the gluon propagator and T u is the quark-gluon vertex. 

QCD-based DSE calculations of F% m (Q 2 ) exist [13, 14], 
the most systematic of which [14] predicted the measured 
form factor [1]. Our primary goal is to elucidate the 
sensitivity of F% m (Q 2 ) to the pointwise behaviour of the 
interaction between quarks. We therefore describe how 
predictions for pion properties change if quarks interact 
instead through a contact interaction. 

We thus begin by defining 



1 



g D^ip-q) T' 
rtvh 



(4) 



where mg is a gluon mass-scale (such a scale is gener- 
ated dynamically in QCD, with a value ~ 0.5 GeV [15]) 
and proceed by embedding this interaction in a rainbow- 
ladder truncation of the DSEs. This means T u (p, q) = 7„ 
in both Eq. (3) and the construction of K in Eq. (2). 
Rainbow-ladder is the leading-order in a nonperturba- 
tive, symmetry-preserving truncation [6, 7]. It is known 
and understood to be accurate for pseudoscalar mesons 
[8, 16] and guarantees current conservation [17]. 
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TABLE I: Results obtained with (in GeV) m = 0, m G = 0.11 , 
Ai r = 0.24 , A uv = 0.823. They are commensurate with those 
from QCD-based DSE studies [20]. Dimensioned quantities 
are listed in GeV or fm, as appropriate, and k := — (gg) 1 . 
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M K f° f° 




r° r° 

'TV ' 7T 




0.23 4.28 0.69 0.68 


0.40 0.22 0.094 0.11 


0.29 0.41 



Preserving the vector and axial- vector Ward-Takahashi 
identities is crucial when computing properties of the 
pion. The m — axial- vector identity states 

P ti T 5 „(k+,k) = S- 1 (k+)h 5 +i 75 S- 1 (k), (12) 

where (&:+,&:) is the axial- vector vertex, which is de- 
termined by 



With this interaction the gap equation becomes 



S 1 (p) = 17 • p + m + 



4 1 

3 TOg 



d 4 q 



7m S(q) 7 M . (5) 



The integral possesses a quadratic divergence, even in the 
chiral limit. If the divergence is regularised in a Poincare 
covariant manner, then the solution is 

Sip)- 1 =i 7 .p + M, (6) 

where M, momentum-independent, is determined by 

Mf 00 , 1 



M = m 



3ir 2 m 2 G j 



s + M 2 



(7) 



r 5/J (fc+,fc) 



4 1 



(13) 

To achieve this, one must implement a regularisation that 
maintains Eq. (12). To see what this entails, contract 
Eq. (13) with P M and use Eq. (12) within the integrand. 
This yields the following two chiral limit identities: 



M 




8 M 
3 m 2 



d 4 q 
(2tt)4 



1 



1 



q 2 + M 2 



d 4 q 



P-q, 



P 



l + M 2 



M 2 



M 2 



(14) 
(15) 



which must be satisfied after regularisation. Analysing 
the integrands using a Feynman parametrisation, one ar- 
rives at the follow identities for P 2 — = m: 



To proceed it is necessary to be specific about the reg- 
ularisation procedure. We write [18] 



1 



M 2 



dTe -r(s+M>) 



dr ( 



-t(s+M 2 ) 



(8) 



where Tj liUV are, respectively, infrared and ultraviolet reg- 
ulators. A nonzero value of r ir =: 1/A ir implements con- 
finement by ensuring the absence of quark production 
thresholds [19]. Furthermore, since Eq. (4) does not de- 
fine a renormalisable theory, A uv := l/r uv cannot be re- 
moved but instead plays a dynamical role and sets the 
scale of all dimensioned quantities. 
The gap equation can now be written 



M = m 



M 



3ir 2 rriQ 



C(M 2 ;T h ,T m ) 



(9) 



where C/M 2 = T(-1,M 2 t 2 v ) - T{-1,M 2 t 2 ), with 
T(a,y) being the incomplete gamma-function. Results 
obtained in the chiral limit are presented in Table I. 

Using the interaction we've specified, the homogeneous 
BSE for the pseudoscalar meson is (q + = q + P) 



d 4 q 
(2n) 4 



(10) 



With a symmetry-preserving regularisation of the inter- 
action, the Bethc-Salpctcr amplitude cannot depend on 
relative momentum. Hence Eq. (1) can be rewritten 



IV(P) = 75 



iK(P) + jjl ■ PF n (P) 



(11) 



M 




16 M 

l 4 q 



d 4 q 



(2tt)4 [q 2 
M 2 



A/ 2 



(2tt) 4 [q 2 + M 2 f 



(16) 
(17) 



Equation (16) is just the chiral-limit gap equation. 
Hence it requires nothing new of the regularisation 
scheme. On the other hand, Eq. (17) is novel: it states 
that the axial- vector Ward-Takahashi identity is satisfied 
if, and only if, the model is regularised so as to ensure 
there are no quadratic or logarithmic divergences. Un- 
surprisingly, these are the just the circumstances under 
which a shift in integration variables is permitted, an 
operation required in order to prove Eq. (12). 

One can now write the explicit form of Eq. (10): 



(18) 



" ' E„{py 
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K-EE 


K-EF 




' E„(P) ' 




" Zn 2 m 2 G 


K-FE 


ICff 




Ftt(P) _ 



where, with m = = P 2 , anticipating the Goldstonc 
character of the pion, 



Kee =C(M 2 ;^ 2 r ,r u 2 v ), K EF 
2JCfe=C 1 (M 2 ;t 2 ,t 2 v ), K F f 



0, 



-2K 



fe ■ 



(19) 



Here C\{z) = —zC'(z), where we have suppressed the 
dependence on Ti r uv . In order to obtain this form of 
ICff, one must employ Eq. (17). The solution of Eq. (18) 
gives the pion's chiral-limit Bcthe-Salpetcr amplitude: 



El = 0.987, F 4 = 0.160, 
written with unit normalisation. 



(20) 
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The canonical normalisation procedure ensures unit 
residue for the pion bound-state contribution to 
the quark-antiquark scattering matrix, a property of 
r£(P) = £ri(P), where 

d 4 q 



N 2 P u = N c ti 



(27T)' 



dPu. 



In the chiral limit, 



P)Ti(P)S(q). 

(21) 



2 2 

'in ' i 



(22) 



The pion's leptonic decay constant is obtained from 
the canonically normalised amplitude and in the chiral 
limit 



J 1 ^ 2 

f. f 



(23) 



If we have preserved Eq. (12), then, in the neighbour- 
hood of P 2 = 0, the solution of Eq. (13) has the form 



P 



(24) 



and the following generalised Goldberger-Treiman rela- 
tions [9] will hold: 

f°E°=M,2^+F R = l. (25) 

That they do can be verified from Table I, which also 
shows that F^{P), necessarily nonzero in a vector ex- 
change theory, irrespective of the pointwisc behaviour of 
the interaction, has a measurable impact on the value of 
/„■, acting here to reduce it by 15%. 

There is more to be learnt from Eqs. (22), (23). These 
chiral-limit expressions exhibit striking similarities. So, 
replace T% by T\ in Eq. (23) and then take the ratio: 

M 



M 2 



M 
El 



(26) 



The last identity is a tautology: it states that E\ = 
M//°, which we already knew from Eq. (12). This illus- 
trates a general result [13]: in a vector exchange theory 
with DCSB, the canonical normalisation constant for the 
pion Bethe-Salpeter amplitude is equivalent to its lep- 
tonic decay constant when both are evaluated in the chi- 
ral limit and a symmetry preserving regularisation pro- 
cedure is employed. 

With the foundation laid, one can evaluate the elec- 
tromagnetic pion form factor in the generalised impulse 
approximation; i.e., at leading-order in the symmetry- 
preserving truncation scheme [13, 14, 17]. Namely, for 
an incoming pion with momentum P \ = K — Q/2, which 
absorbs a photon with space-like momentum Q, so that 
the outgoing pion has momentum p 2 = K + Q/2, 



2K^r{Q 2 ) 



2N r 



dH 
(2^Y 



-tl'D 



iT c n (- P2 )S(t + P2 ) 



xi lfl S(t+ Pl )tTl( Pl )S(t) 



(27) 
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FIG. 1: Fn(Q 2 ) and the separate contributions introduced in 
Eq. (28). 



The form factor is expressible as a sum; viz. 



FTW 2 ) = KIe + K™ef + K m F F , (28) 

= T* E (Q 2 )Ef + T* F (Q 2 )KK + t ff(Q 2 )K 2 - 



In the chiral limit one finds 
N, 



1 EE 



1 EF 



rpTT 

1 FF 



^ 2 JO 
~2Tl E 



dx- 



1 



<j{x) 



Ci(cr(x)) , 



(29) 



--j— o / dxidx2^T-, -C 2 {uj(xi,X2)), (30) 

4tH Jo lj' ! (x-l,X2) 



N c f 1 x\Q 2 
—j / dx\dx 2 



1 



uj{xi,x 2 ) 



Ci (uj(xi,x 2 )) 



uj(xi,x 2 ) + 2AI 

UJ 2 (X1,X2) 



r'2 



C 2 (w(xi,x 2 )) 



(31) 



where: C 2 (z) — (z 2 /2)C"(z) - we have suppressed the 
dependence on Ti I%uv ; a(x) = M 2 + Q 2 x(l — x); and 
u>(xi,x 2 ) = M 2 + Q 2 x\x 2 (\ — x 2 ). It is straightforward 
to evaluate the remaining integrals numerically. 

In Fig. 1 we plot F° m (Q 2 ) and the three separate con- 
tributions defined in Eq. (28). We highlight two features. 
First, F^ n EF (Q 2 = 0) contributes roughly one-third of 
the pion's unit charge. This could have been anticipated 
from Eq. (22). Second, and perhaps more dramatic, is 
that the interaction in Eq. (4) generates 



cw 2 



constant. 



(32) 



Both results originate in the nonzero value of P T (P), 
which is a straightforward consequence of the symmetry- 
preserving treatment of a vector exchange theory [9]. 
From our perspective, Eq. (32) is not a surprise: with 
a symmetry-preserving regularisation of the interaction 
in Eq. (4), the pion's Bethe-Salpeter amplitude cannot 
depend on the constituent's relative momentum. This is 
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FIG. 2: F^ m (Q 2 ) obtained with the interaction of Eq. (4). 
Solid curve - complete calculation with the full Bethe- 
Salpeter kernel; and dashed curve - result obtained with ICfe 
artificially set to zero in Eq. (18). 

characteristic of a pointlike particle, which must have a 
hard form factor. 

It is nonetheless worthwhile to elucidate the math- 
ematical origin of Eq. (32). As has long been known, 
F° n EE has its maximum value at Q 2 = and decreases 
uniformly with increasing Q 2 : for Q 2 ^> M 2 , F°"j^ E oc 
M 2 IQ 2 . However, the mixed contribution, F^"j^ F , is sig- 
nificant in magnitude and negative at Q 2 = 0. As noted 
above, the sum F^ r % E + F^^p provides a pion with unit 
charge. The symmetry-preserving character of our regu- 
larisation scheme has ensured this, as can be seen through 
a comparison of this sum with Eq. (22). Very impor- 
tantly, F^pp increases with Q 2 . It passes through zero 
at Q 2 m 1.5 GeV 2 and continues to grow, approaching a 
nonzero constant value as Q 2 — > oc. These properties are 
apparent from Eq. (30) and follow from the fact that F£ is 
always multiplied by 757 • Q in Eq. (27). F°™p F is zero at 
Q 2 = and evolves to a negative constant as Q 2 — > 00. It 
is always smaller in magnitude than F° n EF . NB. Ref. [13] 
demonstrated that the presence of pseudovector compo- 
nents alters the asymptotic form of F° m (Q 2 ) by a factor 
of Q 2 cf. the result obtained in their absence. 

A striking feature, evident in Fig. 1 and emphasised 
by Fig. 2, is that F° n EE is only a good approximation 
to the net pion form factor for Q 2 < M 2 . F^" EE and 
F° n EF evolve with equal rapidity - there is no reason for 
this to be otherwise, as they are determined by the same 
mass-scales - but a nonzero constant comes quickly to 
dominate over a form factor that falls swiftly to zero. 

In Fig. 3 we compare the form factor computed from 
Eq. (4) with contemporary experimental data [1-3] and 
a QCD-based DSE prediction [14]. Both the QCD- 
based result and that obtained from the momentum- 
independent interaction yield the same values for the 
pion's static properties. However, for Q 2 > the form 
factor computed using ~ 1/fc 2 vector boson exchange is 



FIG. 3: Solid curve - Q 2 F 7r (Q 2 ) obtained with Eq. (4). 
Dashed curve - DSE prediction [14], which employed a 
momentum-dependent renormalisation-group-improved gluon 
exchange interaction. For Q 2 > 0.17 GeV 2 « M 2 , marked 
by the vertical dotted line, the contact interaction result for 
Ftt(Q 2 ) differs from that in Ref. [14] by more than 20%. The 
data are from Refs. [1-3]. 

immediately distinguishable empirically from that pro- 
duced by a momentum-independent interaction. 

Indeed, the figure shows that for F° m , existing exper- 
iments can already distinguish between different possi- 
bilities for the quark-quark interaction. Treated in pre- 
cisely the same manner; i.e., as providing the kernel in 
a symmetry-preserving rainbow-ladder truncation of the 
DSEs, the contact interaction yields a form factor that 
is patently in conflict with the data, whereas the result 
obtained with a momentum-dependent renormalisation- 
group-improved one-gluon exchange interaction agrees 
very well. A dressed-quark propagator and momentum- 
dependent Bethe-Salpeter amplitude based on this in- 
teraction also produces a form factor whose asymptotic 
power-law behaviour agrees with that of perturbative 
QCD [13], unlike the contact interaction. 

In Fig. 4 we illustrate two effects. Firstly, that of a 
nonzero current-quark mass. Namely, with m — 8 MeV: 
m w = 0.14 , M = 0.41 , f„ = 0.0945 (in GeV); and 

El(m n ) = 0.987, ^(ra,) = 0.163. (33) 

The resulting electromagnetic form factor is depicted in 
the figure. Plainly, a nonzero and physical light-quark 
current-mass has no impact on our discussion. 

We also show the impact on F° m of a reduced strength 
for the pion's pseudovector component. For this calcu- 
lation, we suppressed F^(P) by 50%, whilst maintaining 
unit normalisation; viz., we defined 

El = 0.997 , Fl = 0.080 (34) 

cf. Eq. (20), and then recalculated the form factor with 
the pion Bethe-Salpeter amplitude obtained therefrom. 
Notably, this marked suppression of the pseudovector 
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FIG. 4: Dash-dot curve - F^ ul (Q 2 ) computed with = 
0.14 GeV, yielding tv = 0.290 fm; solid curve - chiral limit 
result, rv = 0.292; dotted curve - chiral-limit result obtained 
by suppressing the strength of the pion's pseudovector com- 
ponent by 50% [see Eq.(34)], r n = 0.36 fm; and dashed curve 
- DSE prediction of Ref. [14], tv = 0.67 fm. The data are 
from Refs. [1-3]; and, experimentally, r n = 0.672 ± 0.008 fm 
[21]. 

strength has no material impact on our discussion, even 
though it reduces to just 0.6% the pseudovector compo- 
nent's contribution to the quadratic norm. We emphasise 
that Eq. (34) is not a solution of the pion's BSE. Instead, 
it is an artificial construction employed solely for the pur- 
pose of illustrating the robust nature of our conclusions. 

We explored the consequences for pion observables of 
a momentum-independent vector-exchange interaction, 
regularised in a symmetry-preserving manner. With this 
implementation of the interaction, the results are di- 



rectly comparable with experiment, computations based 
on well-defined and systematically-improvable trunca- 
tions of QCD's DSEs [14], and perturbative QCD. In a 
vector exchange theory, independent of the pointwise be- 
haviour of the interaction, the pion possesses components 
of pseudovector origin, which play a material role in the 
evolution of the pion's form factor away from Q 2 = and, 
indeed, completely determine the form factor's power-law 
behaviour at large Q 2 [13]. 

We find that the contact interaction, whilst capable of 
describing pion static properties well, Table I, generates 
a form factor whose evolution with Q 2 deviates markedly 
from experiment for Q 2 > 0.17 GeV 2 « M 2 and produces 
asymptotic power- law behaviour, Eq. (32), in serious con- 
flict with perturbative-QCD [10-12]. 

The contact interaction produces a momentum- 
independent dresscd-quark mass function, in contrast to 
QCD-bascd DSE studies [19, 22] and lattice-QCD [23]. 
This is fundamentally the origin of the marked discrep- 
ancy between the form factor it produces and extant 
experiment. Hence our study highlights that form fac- 
tor observables, measured at an upgraded Jefferson lab- 
oratory, e.g., are capable of mapping the running of the 
dressed-quark mass function. We are currently working 
to establish the signals of the running mass in baryon 
clastic and transition form factors. 
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